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ABSTRACT 

We present a Bayesian seismology inversion technique for propagating magnetohydrodynamic (MHD) trans- 
verse waves observed in coronal waveguides. The technique uses theoretical predictions for the spatial damping 
of propagating kink waves in transversely inhomogeneous coronal waveguides. It combines wave amplitude 
damping length scales along the waveguide with theoretical results for resonantly damped propagating kink 
waves to infer the plasma density variation across the oscillating structures. Provided the spatial dependence 
of the velocity amplitude along the propagation direction is measured and the existence of two different damp- 
ing regimes is identified, the technique would enable us to fully constrain the transverse density structuring, 
providing estimates for the density contrast and its transverse inhomogeneity length scale. 

Subject headings: magnetohydrodynamics (MHD) — methods: statistical — Sun: corona — Sun: oscillations 



1. INTRODUCTION 

Observations show that propagating magnetohydrodynamic 
(MHD) waves are present in the solar atmosphere. Small am- 
plitude propagating transverse MHD waves have received par- 
ticular attention, because of their potential for energy storage, 
transfer, and deposition in the context of wave-based plasma 
heating processes. A key problem yet to be solved is the quan- 
tification of the role of waves in coronal heating. The solu- 
tion requires a reliable knowledge of the physical properties 
in magnetic and plasma structures. Because physical proper- 
ties, such as the magnetic field strength, the plasma density, 
and their field-aligned and cross-field structuring cannot be 
directly measured, seismology of MHD waves offers an alter- 
native method to probe coronal plasmas. 

The presence of ubiquito us coronal t ransv erse waves 
has been reported by Tomczv k et al.l (J2007) and by 
lTomczyk& Mcintosh (2009) in observations taken with the 
Coronal Multi-Channel Polarimeter (CoMP). The distur- 
bances have amplitudes of the order of 0.3 km s ' and prop- 
agate at speeds of about 0.6 megameters per second. They 
are interpreted as transverse MHD waves of Alfvenic charac- 
ter The measured discrepancy betwee n inward and outward 
powe r associated to the disturbances (iTomczvk & Mclntoshl 
120091) is thought to be an indication of in situ wave damp- 
ing. Resonant absorption of kink waves has been proposed 
to explain this damping bv'Pascoe et alj (120101) . This process, 
studied in the context of standing ki nk waves in coronal loops 
(see e.g. IG oossens et al. 1992; Ruderman & Roberts 20021 
iGoossense t al. 2002; Goossens et al. 2006) and in pro mi 
nence fine structures (Arregui et al. 2008; S oler et 



n promi- 
al.l 120091) 



occurs because of transverse inhomogeneity of the plasma 
across the waveguides. Resonant absor ption predicts a selec - 
tive damping as a function of frequency (iTerradas et al.l20foh . 
The good agreement between the observationally measured 
and the theoretically predicted outward/inward power ratio 
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as a function of wave frequency strongly suppo rts the idea 
that r esonant absorption operates on these waves (IVerth et alJ 

lIoToh . 

The use of wave damping properties to perform seis- 
mology inver s ions w ith transverse wave s was proposed by 
lArregui et d] (120071) and IGoossens et al.l (|2008") in the con- 
text of standing k i nk wa ves in coronal loops, and by 
lArregui & Ballesten (120111) i n the context of prominence 
fine s t ructures (see reviews by IGoossens! 120081 lArregui et al.l 
120 12t lArreguil 12012). Because the number of unknowns 
is usually larger than that of observables, density mea- 
surements need to be used as additional observational in- 
formation to fully constrain the parameters of interest 
(lArregui & AsensioRamos 20 1 1 ). 

The mechanism of resonant absorption does not make a dis- 
tinction between standing and propagating MHD kink waves. 
The dynamics correspon ds to a surface Alfven wave d amped 
by resonant absorption ( Goossens et al.l l2009l l2012al) . For 
standing kink waves, the observational consequence is the 
attenuation of the amplitude in time. For propagating kink 
w aves, the attenuation of wave amplitude in space. 

IGoossens et al.l (l2012bl) have presented a seismology inver- 
sion scheme for propagating MHD waves damped by resonant 
absorption. In their analysis, a solution curve that relates the 
Alfven speed, the density contrast and the transverse inhomo- 
geneity length scale is obtained. The Alfven speed is con- 
strained to a narrow range, but the two parameters that define 
th e cross-field d ensity structuring cannot be inferred. 

iPascoe et alj {2012) have pointed out that a more general 
damping profile consisting of a Gaussian and an exponen- 
tial profile better reproduces the spatial decay of the ampli- 
tude for propagating kink waves. A Gaussian seems to bet- 
ter reproduce the amplitude behavior initially, while the ex- 
ponential profile properly describes the dynamics after a few 
wavelengths. Their numerical solutions point to a transition 
between the two damping regimes. Currently, only indirect 
observational evidence about in situ damping of propagating 
transverse waves is available, with no measurement of the 
damping length-scales. The confirmation on the existence of 
two damping regimes would make accessible information on 
two damping length scales and the height at which the transi- 
tion between the two regimes occurs. 
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In this Letter, we demonstrate how this information can be 
used to fully constrain the cross-field density structuring in 
coronal waveguides. 

2. SPATIAL DAMPING OF PROPAGATING KINK WAVES 

The theory for the spatial damping of propagating kink 
waves due to resona nt absorption has been developed by 
a number of studies (Terradas et al."2010t IVerth etaP 120101 
[Soler et al. 201 la c.b; Hood et al. 2013). Numerical simula- 
tions have confirmed the obtained damping properties, by an- 
alyzing spatial and temporal properties of the mode coupling 
process in coro nal loops and arbit r ary inhomogeneou s coro- 
nal sti-uctures ( Pascoe et all l20Toi |20TTI l20T2ll20T3l) . The 
classic theoretical model assumes that transverse kink waves 
are guided by plasma structures. The waveguides are con- 
sidered to be density tubes in the zero plasma-/? approxima- 
tion with a uniform magnetic field directed along the axis of 
a straight cylindrically symmetric structure. They have a uni- 
form internal density, p,, and a uniform external density, pe, 
with Pi > Pf , connected in the transverse direction by a non- 
uniform region of thickness IjR, with R the tube radius. Be- 
cause of the non-uniformity, resonant absorption produces the 
decay of the wave amplitude in space as the wave propagates. 

For propagating waves with a real frequency co, an analyti- 
cal expression for the wavelength in the thin tube approxima- 
tion is (see.Terrad as et al. :, 201 0) 



A = V2vA,ir 
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Here va.i is the internal Alfven velocity, T the period, and 
( - pi/pe the density contrast. This solution is valid provided 
(L)R/v\,i « 1. It expresses a relation between two observable 
quantities (wavelength and period) and two quantities to be 
determined (Alfven velocity and density contrast). 

Because of resonant absorption, spatial damping occurs, 
and the transverse velocity amplitude decays with an expo- 
nential profile of the form exp(-z/Ld). Under the thin tube 
and thin boundary (l/R « 1) approximation, an expression 
for the damping length, Ld, as a function of the relevant phys- 
ical parameters can b e obtained. In units of the wavelength 
this expression is (see lTerradas et al.]|2010b 
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The first factor is due to the assumed linear density profile at 
the non-uniform layer Note that the right hand-side of this 
expression is identical to the one for the damping time over 
the period for standing kink waves. The reason is that reso- 
nant absorption does not make any distinction with respect to 
the standing or propagating character of the wave. 

The exp onential profile obtai ned for standing (e.g., 
iRuderman & Roberts 2002; Gooss ens et al.l 12002) and prop- 
agating (e.g.. iTerradas et al.,i2010.) kink waves describes the 
asymptotic s tate of the damping behavior, i.e. at large times 
or distances. iPascoe et al.l (1201 2*) demonstrated with numeri- 
cal simulations that the initial damping stage can be described 
by a Gaussian profile of the form exp(-z^/L g), wit h Lg the 
Gaussian damping length scale. iHood et al] ( 120131) consid- 
ered the problem analytically and produced an expression for 
the full nonlinear spatial damping profile, which can be ap- 
proximated as Gaussian for low heights and exponential at 




Figure 1. Transverse velocity component as a function of height at the axis 
of the tube for propagating kink waves for a numerical simulation with f = 
1.5 and l/R = 0.4 (case 4 in Table |2)- On the top is the general spatial 
damping profile given by the sohd line. The transition between Gaussian 
and exponential damping is given by the vertical dotted line. On the bottom, 
the general profile is split into its two components; Gaussian (dot-dash) and 
exponential (dashed). 



large heights. Instead, iPascoe et al.l ( 120 13h proposed a gen- 
eral spatial damping profile composed of a Gaussian damp- 
ing profile at low heights and the usual exponential profile at 
large heights. An example of the spatial dependence of the 
velocity amplitude from numerical simulations and the dou- 
ble profile fitting for such a general damping profile is dis- 
played in Figure[T] The accuracy of this approximate damping 
pr ofile was demonstrat ed by the parametric study performed 
bv IPascoe et alj (120131) . This study shows that the Gaussian 
damping length scale can be well described by the expression 
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This equation expresses the Gaussian damping length as a 
function of the same two parameters that determine the ex- 
ponential damping length. This means that the observational 
identification of two damping regimes and the measurement 
of their associated length scales would provide us with addi- 
tional information without the inclusion of new model param- 
eters. The height, h, at which the damping regime c hanges 
from Gaussian to exponential is given by (see .Pascoe et al.l 
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The non-linear (non-exponential) damping rate is an inher- 
ent property of resonantly damped kink oscillations, as de- 
rived a nalytically by Hoodet al. (2 013h for spatial damping 
and by Ruderm an & Terra das (2013) for time damping. The 
previous analyse s by e.g., Ruderman & Roberts (2002) and 
iGoossens et alJ ((2002), only attempt to describe the damping 
in the asymptotic state and so they only obtain the later expo- 
nential damping stage. 

The observational identification of this feature would be 
strong support for the resonant damping mechanism. How- 
ever, no measurement of the axial damping spatial scales is 
available yet. Our analysis shows that, if present. Equa- 
tions (|2]l-(|4li can be used to fully determine the transverse 
density structuring in oscillating coronal waveguides. 

Out of the three equations, only two are independent. In 
our inversion scheme, Eqs. (|3) and (|4) will be used to in- 
fer ^ and l/R from data given by Lg and h. The two alge- 
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braic equations can be directly solved for the two unknowns. 
However, the problem we intend to solve is not exact, since 
real data are noisy and measurements have an associated 
uncertainty. To accommodate a proper propagation of un- 
certainty, we make use of Bayesian analysis, a framework 
that is starting now to be developed in coronal seis mology 
JArregui & Asensio Ramosll201 UlArregui et aLll2013l) . 

3. BAYESIAN INVERSION SCHEME 

To perform our inversion we use Bayes' theorem 
(lBaves&Pricelll763h 



p(0\d) = 



p{d\e)p{e) 

p(d) ' 



(5) 



which gives the solution to the inverse problem in terms of 
the posterior probability distribution, p(6\d), that describes 
how probability is distributed among the possible values of 
the unknown parameter, 9, given the data d. The posterior is 
a combination of what is known about the parameters before 
taking into account the data, the prior p{0), and the likelihood 
function, p{d\0), which tells us how well the model predicts 
the observed data. The denominator, p{d), is the so-called ev- 
idence, which normalises the likelihood. It plays no role in 
parameter inference. 

Once the posterior is known, we calculate how a particular 
parameter is constrained by data computing its marginal pos- 
terior. This is done by performing the following integral of 
the full posterior with respect to the rest of parameters 



p(0i\d) 
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The result provides us with all the information for model pa- 
rameter 0i available in the priors and the data. This method 
also enables us to correctly propagate uncertainties from data 
to inferred parameters. 

We next specify the likelihood function and the priors. 
In what follows we assume the observed data are given by 
d - (Lg, h), where both observed length-scales are normal- 
ized to the wavelength. The unknowns are gathered in the 
vector 0-{^, l/R). Under the assumption that observations are 
corrupted with Gaussian noise and they are statistically inde- 
pendent, the likelihood can be expressed as 



p(d\0) = [Ino-L^o-hj exp j — -i- 



, .. 2cr2 J' 

(7) 
with L'^^iO) and h'y\0) given by Equations Q and ©. Like- 
wise, cr^ and aj^ are the variances associated to the Gaussian 

damping length and the height h, respectively. 

The priors indicate our level of knowledge (ignorance) be- 
fore considering the observed data. We have adopted uniform 
prior distributions for both unknowns over given ranges, so 
that we can write 



p(Od = 
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for 6^'" < < 0™ 
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and zero otherwise. For the minimum and maximum values 
the intervals ^ e [1, 20] and l/R e [0, 2] have been taken. This 



choice of priors expresses our belief that the unknown param- 
eters are constrained to those ranges, with all combinations 
being equally probable. 

Computations were repeated for different uniform prior 
ranges for f, with ^min = 1 and ^max - 5, 10,50, and for 
l/R, with {l/R)mm = and (l/R)^^^ = 1, 1.5,2. We also used 
strongly informative priors (Gaussians in ^ and l/R centred 
far away the synthetic parameter values). For strong prior 
information on l/R, convergence to the correct posterior is 
achieved through repeated application of Bayes' theorem, us- 
ing the posterior of one iteration as a prior for the next one 
and multiplying again by the likelihood. This means that a 
few consistent observations of the data lead to a converged 
posterior, independently of the prior Our sensitivity analy- 
ses lead us to conclude that posteriors are dominated by the 
information contained in the data, that overwhelms the prior 
information. 

Posteriors are evaluated for different combinations of pa- 
rameters using Bayes' theorem (|5) and inferred distributions 
are obtained from the computation of the marginal poste- 
riors using Eq. (|6]l. To solve the ID integrals, numeri- 
cal quadratures were employed, using an adaptive Gauss- 
Kronrod quadrature. All computations were checked by com- 
pa rison to results obtained with the MCMC code employed 
bv lArregui & Asensio Ramosl(l201 II) . 

4. INVERSION RESULTS 

We first evaluated the performance of our inversion scheme 
by making the inference under controlled conditions. We gen- 
erated predictions for the length scales Lg and h for different 
combinations of the equilibrium parameters, ( = 1.5,2,3,4 
and l/R = 0.05,0.15,0.2,0.4, using Equations Q and ©. 
Those synthetic data where treated as observed data in the 
Bayesian inversion. A 10% uncertainty on the data was con- 
sidered and the posterior distributions for ^ and l/R were com- 
puted, using the likelihood function (JTl and the uniform pri- 
ors. Once the posteriors were known, the median and the vari- 
ances associated to the 68% confidence level were calculated. 
Table[T]displays the inversion results for some parameter com- 
binations. In all the cases, the inversion method correctly re- 
covers the values for the physical parameters. The larger the 
density contrast, the shorter the two length-scales Lg and h 
are. This increases the errors in the inferred density contrast, 
while errors in l/R are not affected that much. For the combi- 
nations with the largest ^ = 10 and l/R = 0.5, 1, 1.5, L„ and 
h are comparable to the wavelength. This would make very 
problematic the observational identification of the two damp- 
ing regimes. 

Then, simulations of transverse kink wave propagation in 
a m agnetic flux tube w ere performed using a numerical code 
(see lPascoe et al.ll2013 , for details). A Lax- Wendroff code is 
used to solve the linear MHD equations in cylindrical coordi- 
nates. The lower boundary is driven harmonically with veloc- 
ity perturbations corresponding to the loop footprint moving 
back and forth about its equilibrium. The simulation ends af- 
ter 10 periods of oscillation and the spatial damping profile is 
investigated by considering the radial velocity component, v^ 
as a function of z at the centre of the loop (Figure [T). From 
the behavior of the amplitude of the excited kink waves at dif- 
ferent heights the damping profile was fitted and values for Lg 
and h obtained. Using those fitted values as data, we repeated 
the inversion procedure. For the sake of comparison, parame- 
ter spaces that overlap with those in Table [T] were considered. 

Figure|2]displays an example of the marginal posterior dis- 
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Figure 2. One-dimensional marginalized posterior distributions for tlie density contrast (a) and the transverse inhomogeneity length scale (b) corresponding to the 



inversion of a spatially damped transverse oscillation with Lg/A = 5.0± 0. 1 and h/A -- 
posterior distribution. The light and dark grey shaded regions cover the 95% and 68"; 



4.9 ±0.1 (uncertainties of 10% have been used), (c) Joint two-dimensional 
credible regions. The symbol indicates the estimate. 



Table 1 

Inversion of Synthetic Data Using the Analytical Forward Model 



Syntht 


tic Parameters 
l/R 


Synthetic Data 
L,/A h/A 


Inversion Results 


i 


f 


l/R 


1.5 


0.05 


14.2 


5.0 


1 ^,+0.08 


05+0''- 


1.5 
1.5 


0.15 
0.2 


8.2 
7.1 


5.0 
5.0 


1.50+""'' 
1.51+§|^ 


0.1 e+o"^ 
o.2i+o:| 


1.5 


0.4 


5.0 


5.0 


1 50+°-''' 
'•'"-0.05 


44+0:13 


3 


0.05 


5.7 


2.0 


3.11+«-5' 


0.05+002 


3 


0.15 


3.3 


2.0 


3-09 oi 


0.15+™= 


3 


0.2 


2.9 


2.0 




19+™* 


3 


0.4 


2.0 


2.0 


3.10+°i? 


0.42+0-15 


4 


0.05 


4.8 


1.7 


^■-"-0.79 


o'i+oM 


4 


0.15 


2.7 


1.7 


4.39+ig 


0.15+o°5 


4 


0.2 


2.4 


1.7 


4.38+I-I 


iQ+o:o8 


4 


0.4 


1.7 


1.7 




n ao+0.14 


10 


0.5 


1.1 


1.2 


51+01* 
"■■"-0.11 


10 


1.0 


0.8 


1.2 


11.55!^-^? 


1 02+0-29 
^■^-^-0.28 


10 


1.5 


0.6 


1.2 


1 7 70+4.32 



tributions and the joint probability distribution for ( and l/R. 
For both parameters, well defined probability distributions are 
obtained. For each parameter, the median of the marginal pos- 
terior and errors given at the 68% credible interval are used 
to compute the estimates given in Table |2l This Table shows 
the values for the physical parameters used in the simulations. 
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Inversion of Numerical Data From Simulations 




Simula 


tion Parameters 
l/R 


Fitted Data 

L,/A h/A 


Inversion Results 


i 


f 


l/R 


1.5 


0.05 


11.5 


3.8 


1 70+0.12 


Q Q^+0.02 


1.5 


0.15 


7.9 


4.6 


1.56+008 


0.15+005 


1.5 


0.2 


7.0 


4.8 


1-53 06 


0.21+0:0] 


1.5 


0.4 


5.0 


4.9 


1.52+0:07 


0.39+0™ 


3 


0.05 


5.5 


2.1 


2.88+o:| 


0.06+o|2 


3 
3 


0.15 
0.2 


3.5 
3.1 


2.2 
2.2 


274+o:w 

2 74+8:t? 

Z./4„3Q 


0.16+o:oB 
A2i+S:8f 

U.Z1_Q05 


3 


0.4 


2.1 


2.0 


3.09+0" 


0.38+0}] 


4 


0.05 


4.9 


1.7 


4-17tf 


0.05+0:0^ 


4 
4 


0.15 
0.2 


3.1 

2.7 


1.9 
1.9 


3 J 9+0:64 
3.33+o:| 


0.16+o:o6 
0.21+|§] 


4 


0.4 


2.3 


2.2 


2 73+033 


0.38^ 



the fitted length scales, and the inferred physical parameters. 
Numerical and analytical forward models give similar results. 
This issue is discussed in detail by Pascoe et al. (2013) (see 
their figures 8, 9, and 10). Our Bayesian inference technique 
properly returns the physical parameters of interest. This 
means that the algebraic expressions for the forward problem 
are accurate enough, since their use correctly returns the val- 
ues for the physical parameters that were actually used in the 
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numerical setup. As with synthetic data in Table [T] large den- 
sity contrast values tend to produce larger errors in their de- 
termination by inversion. The main problem lies in obtaining 
the parameters Lg, Lj, and h from the data, and specifically in 
determining h accurately, which determines the accuracy of 
the density estimate. 

The general spatial damping profil e remains an accurat e de- 
scription of the damping behavior (iPascoe et al.l 12013b . for 
large density contrasts. The weak link is the errors in the least 
squares fit of the damping profile to the data. The errors are 
larger for higher density contrasts due to the Gaussian stage 
being very short, and the results are more sensitive to the er- 
rors in h. The lesson is that large contrasts are a challenge 
from the o bservational point o f view, if the t heor etical pre- 
dictions by iHood et all (120 13h : iPascoe et all 12013) and the 
inversion technique presented here are to be employed. 

5. CONCLUSION 

The cross-field density structuring of magnetic waveguides 
cannot be currently estimated using seismology techniques. 
Yet such an understanding is crucial to assess the possible role 
of MHD waves in plasma heating processes. Recent theoreti- 
cal results for the spatial damping of propagating MHD kink 
waves predict the existence of two different damping regimes 
for the spatial dependence of the velocity amplitude along the 
waveguides. In this Letter, we have explained how the ob- 
servational identification of those regimes and the measure- 
ment of the associated damping length scales can be used to 
fully constrain the cross-field density structuring of magnetic 
waveguides. 

Such measurements are currently unavailable and observa- 
tional efforts towards the observational confirmation of the 
theoretical prediction and the measurement of the associated 
spatial scales should be devised. They would not only pro- 
vide us with hitherto undetermined information on density, 
but would strongly endorse resonant absorption as damp- 
ing/heating mechanism, since the presence of two damping 
regimes is an inherent property of this mechanism. 

With the use of Bayesian inference, such measure- 
ments would ensure the inverse problem is consistently 
solved taking into account all the relevant information 
and with correct propagation of uncertainty. The range 
of parameters used in our inversion experiments covers 
the reasonable range of expected values u sed in previous 
studies (Goossens et al. 2002; Ruderman & RobertsI 120021; 
lArregui et al.ll2007t ;Goossens et al.l 120 12 b). so the measure- 
ment and inversion of observations seem feasible. 

The application to real observations would require high 
quality data, so that a reliable fitting can be performed to 
obtain the damping length scales and to be able to clearly 
discriminate between Gaussian and exponential damping 
regimes. As pointed out by Pascoe et al. (2013), the Gaussian 
damping stage only applies for one wavelength, in the limit 
of density contrast tending to infinity. Fitting to such a small 
signal can lead to significant error bars. The method of fitting 
both damping regimes is applicable for low density contrast 
waveguides. Our application is limited to infe rence for a sin- 
gle oscillating magnetic tube. Observations by lTomczyk et al.l 
O007) and Tomczyk & Mcintosh (2009) show that wave dy- 
namics is ubiquitous over large coronal regions. Our method 
could be extended in order to apply the inversion technique 
to such extended regions, thus obtaining overall information 
about the cross-field density structuring of the corona where 
these waves propagate. 



Recent studies have considered effects on the wave ampli- 
tude of propagating transverse kink waves from the consider- 
ation of a background field-aligned flow (ISoIer et al. 1201 Ibl) 
and of longitudinal density stratification (ISoler et al.l 1201 id) . 
These two effects affect the wave amplitude and might par- 
tially or totally compensate the damping by resonant absorp- 
tion, thus producing significant differences in the inferred pa- 
rameters obtained in this study. Future seismology schemes 
should include these ingredients in their inversion process. 
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